Clusters, columns, and lamellae — minimum energy configurations in core softened 

potentials 



00 

o 
o 

<N 

X> 
0) 

00 
<N 



O 

a 

o 
o 



> 

ON 

On 

^t" 

<N 
O 
00 

o 



Gemot J. Pauschenwein 1 and Gerhard Kahl 1 

1 Center for Computational Materials Science and Institut fur Theoretische Physik, 
Technische Universitat Wien, Wiedner Hauptstrafie 8-10, A-1040 Wien, Austria 

(Dated: February 28, 2008) 

We give evidence that particles interacting via the simple, radially symmetric square-shoulder po- 
tential can self-organize in highly complex, low-symmetry lattices, forming thereby clusters, columns, 
or lamellae; only at high pressure compact, high-symmetry structures are observed. Our search for 
these ordered equilibrium structures is based on ideas of genetic algorithms, a strategy that is char- 
acterized by a high success rate. A simple mean- field type consideration complements these findings 
and locates in a semi- quantitative way the cross-over between the competing structures. 
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"One of the continuing scandals in the physical sci- 
ences is that it remains in general impossible to predict 
the structure of even the simplest crystalline solids from 
a knowledge of their chemical composition" [l|. Even 
nowadays, twenty years later, this statement is still valid 
and it applies equally well to the case that the physical 
properties of the system, e.g., in terms of its inter-particle 
potential, are known. Indeed, in many problems of hard 
and soft condensed matter theory a powerful tool that 
is able to predict the ordered equilibrium structures of a 
system in a reliable way is badly missing. This applies in 
particular to soft matter systems where particles are able 
to self-organize in a broad variety of unexpected and of- 
ten very exotic ordered structures: micellar and inverse 
micellar structures 0, 0], spirals chains and layers 
0,01, and cluster phases 0,0] are a few examples. 

During the past decades several strategies have been 
proposed to find the energetically most favorable parti- 
cle arrangements of a system. Apart from conventional 
approaches that rely on intuition, experience, or plau- 
sible arguments when selecting candidates for ordered 
equilibrium structures, there are more sophisticated ap- 
proaches such as simulated annealing [8[, basin hopping 
or meta-dynamics [10[. However, all these strategies 
are affected by different sorts of deficiencies which can 
significantly reduce their success rates. 

In recent years convincing evidence has been given that 
search strategies based on ideas of genetic algorithms 
(GAs) are able to provide a significant break-through to 
solve this problem. Generally speaking, GAs are strate- 
gies that use key ideas of evolutionary processes, such as 
survival of the fittest, recombination, or mutation, to find 
optimal solutions for a problem [ll| . The wide spectrum 
of obviously successful applications in different fields of 
condensed matter physics unambiguously demonstrates 
their flexibility, reliability, and efficiency: among these 
are laser pulse control [12], protein folding [13[, or cluster 
formation [14]. In contrast, attempts to apply GAs in the 
search for ordered equilibrium structures in condensed 
matter theory were realized considerably later. While 
the first application probably dates back to 1999 [l5[ 



their widespread use in hard matter theory was pioneered 
by Oganov and co-workers [l6[ only recently, where they 
have become meanwhile a standard tool: a wide spec- 
trum of successful applications ranging from geophysical 
to technologically relevant problems give evidence of the 
power and the flexibility of this approach (for an overview 
see [16[). I n so ft condensed matter theory the usage of 
these search strategies is still in its infancy. First appli- 
cations to find minimum energy configurations (MECs) 
of soft systems have, nevertheless, unambiguously doc- 
umented the power of the algorithm: successful exam- 
ples are the identification of exotic lattice structures and 
cluster phases for particular soft systems 0, Q3? or °f 
complex, ordered arrangements of monolayers of binary 
dipolar mixtures [18[ . All these investigations mentioned 
above give evidence that GA-based search strategies have 
an extremely high success rate. 

In this contribution we consider a simple soft model 
system, i.e., a square-shoulder system, and show that the 
GA is an efficient and reliable tool to identify even highly 
complex MECs in soft matter systems. We discover an 
overwhelming and undoubtedly unexpected wealth of or- 
dered MECs, which comprise clusters, columns, lamellae, 
as well as compact structures. The simple shape of the 
potential allows an easy geometric interpretation of the 
system's strategies to form MECs in terms of overlapping 
shoulders, a nice feature that is not available in systems 
with continuous potentials. We find evidence that the 
success rate of our algorithm must be close to 100%. 

The square-shoulder system is the simplest represen- 
tative in the class of the so-called core-softened poten- 
tials (for references see (l9[). Despite their innocently 
looking potentials these systems are characterized by a 
host of surprising properties, such as water-like anoma- 
lies [19|, [20| or a rich wealth in the occurring structures, 
investigated in detail in two-dimensional systems @, [l8[ . 
The potential 3>(r) of the square-shoulder system, 
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consists of an impenetrable core of diameter a with an 
adjacent, repulsive, step-shaped shoulder of height e and 
width Act. For the aims of the present contribution it rep- 
resents the 'quintessential' test system [21]. But we also 
point out that the square-shoulder system is not only 
of purely academic interest: it represents a reasonable 
model system for colloidal particles with a core-corona 
architecture, as has been given evidence for in [22[. Fur- 
ther, we introduce the number-density g = N/V, N and 
V being the number of particles and the volume, respec- 
tively. 

The MECs for this system have been identified via a 
GA-based search strategy in the NPT ensemble. For 
details of the encoding of the individuals (= lattices) 
and of the reproduction mechanism we refer to (23[. 
The quality of an individual X is measured via the fit- 
ness function /(X), for which we haven chosen f(T) — 
exp{ — [G(T) — Go] /Go}- Since our investigations are re- 
stricted to the case T = 0, the Gibbs free energy, G(Z), 
reduces at a given pressure P to G = U + PV, U being 
the lattice sum; Go is the Gibbs free energy of a reference 
structure. Significant extensions of the search strategy 
were required due to the impenetrable core: as a con- 
sequence of the highly stochastic character of GAs, the 
algorithm will propose with a high probability lattices 
where the cores of the particles overlap: such configu- 
rations are unphysical and thus useless. A more quan- 
titative investigation reveals that the physically relevant 
lattices (characterized by non-overlapping cores of the 
particles) populate only a highly porous subset of the 
search space |24j . 

To overcome this problem a suitable modification is ur- 
gently required, which guarantees that unphysical config- 
urations are excluded a priori so that only lattices with 
non-overlapping cores are created. Such a strategy has 
been developed and will be outlined briefly. We start 
for simplicity with a simple lattice, described via a set 
of linearly independent vectors ai, a2, and a3. In or- 
der to satisfy our expectations, they have to meet several 
requirements. First, the vectors are chosen such that 
I < |a_2 1 < |a.3 1 , where |ai| represents the shortest pos- 
sible distance between lattice sites. Then, dL2 is selected 
so that |a.2 1 is either equal to |ai| or represents the sec- 
ond smallest distance encountered between two lattice 
points. Finally, the third vector, a3, is chosen so that 
|a_3 1 satisfies a similar relation with respect to |a.2 1 - Of 
course, the choice of the {a^} is not unique. If we are 
able to construct a lattice so that a priori |ai| > <r, then 
it is guaranteed by construction that the particles will 
not overlap. We have succeeded in developing a formal- 
ism that is able to create vectors {a^}, that satisfy the 
above requirements. These constraints lead to inequal- 
ities between the Cartesian components of the vectors; 
these lengthy relations along with a detailed explanation 
of the algorithm will be described elsewhere [25[ . In non- 
simple lattices overlap can also be caused by the basis 
particles. Let us assume that the underlying simple lat- 
tice satisfies above requirements. Then we calculate all 



the distances between the particles of this cell and all 
particles located in the 26 neighboring cells. Let Iq be 
the smallest among these distances; if Iq < |ai|, then the 
vectors {a^} are scaled by a factor |ai|//o, which guar- 
antees that in the entire lattice no particle overlap will 
occur. 

Thermodynamic properties will be presented in stan- 
dard reduced units: g* = gcr 3 , P* = P<r 3 /e, U* = 
U/(Ne), and G* = G/(Ne) = U* + P*/q*. The simple 
shape of the potential allows us to simplify the search 
considerably: since for a fixed particle configuration U* 
(counting the number of overlapping coronas) is con- 
stant, G* = U* + P* I g* is a straight line of slope 1/g* 
in the (G*, P*)-diagram. The limiting MECs at low 
and high pressure are easily identified as close-packed 
spheres, either with diameter Act (corresponding to a 
slope 1/^min = A 3 />/2) or with diameter a (correspond- 
ing to a slope 1/^max = VV^) - cf. Fig. [TJ Any other 
MEC occurring in this system is characterized by a line 
of slope 1/g*, satisfying l/g^ in > 1/g* > 1/^ax- Tnus 
q* _ Q*(p*^ jg a sequence of straight lines. In a first step 
we determine in the (G*, P*)-diagram the intersection 
point of the two limiting straight lines mentioned above 
and launch for this P*-value a GA-search. This leads to 
a new MEC, characterized by a lower G*- value, a density 
g* , and thus by a new straight line of slope 1/g*. Inter- 
secting this line with the previous lines and iterating this 
strategy, the entire pressure-range can be investigated in 
an extremely efficient way. From the economic point of 
view this systematic procedure is very attractive; but it 
has two additional advantages: (i) there is no risk to 'for- 
get' MECs which easily occurs when working on a finite 
P*-grid; and (ii) it avoids a characteristic drawback of 
GAs that tend to fail in the vicinity of state points that 
are characterized by competing structures; in the present 
approach these transition points are determined exactly 
via the intersection of two straight lines. 

At each of the intersection points, 1000 to 3000 inde- 
pendent GA-runs with 700 generations, each consisting of 
500 individuals, have been performed. Up to twelve par- 
ticles per unit cell were considered, a number which of- 
fers the system sufficient possibilities to form even highly 
complex structures. Bearing all this in mind, we are con- 
fident that the sequence of MECs which we present in 
the following are complete. 

We have considered three values of A, corresponding 
to small (A = 1.5), intermediate (A = 4.5), and long 
(A = 10) shoulder width. For the first case, G* and U* , 
as functions of P*, are presented in Fig. [TJ G* is, as 
mentioned above, a sequence of straight lines, while the 
energy levels of U* are rational numbers, given by the 
number of overlapping coronas per particle. For A = 1.5, 
seven MECs can be identified (cf. acronyms in Fig. [1]). 
At low pressure a columnar structure is identified, while 
for larger P*-values, the relatively short corona width 
allows only for compact structures. 

For A = 4.5, the considerably broader corona offers 
the system a large variety of strategies to form MECs. 
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FIG. 1: (color online). G* and U* as functions of P* for 
the square-shoulder system, A = 1.5. Shaded areas represent 
contributions to G* as labeled. Broken lines indicate low- and 
high-pressure limiting cases of MECs (see text). 



Their total number amounts to 33 (cf. Fig. [2j) which can 
be grouped with increasing pressure into four structural 
archetypes: cluster phases, columnar structures, lamellar 
phases, and, finally, compact structures (see also Fig. [4]). 
At low pressure the system forms cluster crystals, i.e., pe- 
riodic structures where the lattice points are populated 
by clusters of particles. A closer analysis reveals a strong 
interplay between the cluster shape and the cell geome- 
try in the sense that overlaps of coronas of neighboring 
clusters are avoided. Inside the clusters, which contain 
up to eight particles, the cores are in direct contact. For 
demonstration we have depicted in Fig. [31(a) a typical 
cluster crystal: the equilateral triangle-shaped clusters 
populate the lattice positions of a single face centered 
monoclinic structure. As the pressure is increased the 
system develops a new strategy to minimize G by forming 
columnar structures. Particles aligned in close or direct 
contact form single or double stranded columns, which 
are arranged periodically in space. The inter-columnar 
spacing is imposed via the width of the corona. An ex- 
ample for the columnar phase is given in Fig. [3^b) . As 
the pressure increases further, the columns are squeezed 
together in one direction, forming thereby lamellae; the 
transition scenario is depicted in Fig. 03c). The lamellar 
structures are realized in a large variety of morphologies. 
In Fig. 02(d) we show an example for a lamellar phase: 
here, two hexagonally close-packed, planar layers are in 
direct contact and form double layers. Finally, under the 
influence of the increasing pressure, the lamellae merge, 
creating in this way typical compact structures. 

For A = 10, the variety of structures is even richer, in 
total 48 MECs have been identified. The clusters contain 
more particles and the columns are more complex in their 
morphology. For demonstration we depict two of these 
structures in Figs. 03(e) and (f). However, the sequence 
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FIG. 2: (color online). G* and U* as functions of P* on a 
double-logarithmic scale for the square-shoulder system, A = 
4.5. Shaded areas represent contributions to G* as labeled. 
Arrows cf. Fig. [3] 



of structural archetypes that has already been identified 
for A = 4.5 is maintained. 

The above mentioned transition from clusters to com- 
pact structures can be understood on a simple, semi- 
quantitative level. Generalizing the ideas proposed by 
Glaser et al. [3] for two dimensional systems, we only 
consider aggregates (i.e., clusters, columns, and lamellae) 
instead of the individual particles. Assuming an ideal- 
ized shape for these aggregates (spheres, straight lines, 
and planes) the inter- and intra- aggregate energy can be 
calculated. The respective expressions are in some cases 
rather complex and will be presented elsewhere [25[. G* 
is then characterized by g*, the distance between two 
aggregates, and their spatial extent. Retaining parame- 
ters up to first order and minimizing G* with respect to 
these quantities we obtain the thermodynamic properties 
of the respective phase. In particular, (U* + l/2)/A 3 as a 
function of P*/A 3 turns out to be A-independent and has 
been plotted, along with the respective results obtained 
via the GA, for all A- values in Fig. [H As expected, agree- 
ment improves considerably with increasing A. A detailed 
analysis of the MECs identified by the GA (in particu- 
lar for the larger A's) reveals, that the MECs populate, 
according to their aggregate- type, nearly exclusively the 
respective pressure-range. 

Finally, a nice by-product should be mentioned: the 
well-defined range of the shoulder represents a highly sen- 
sitive antenna to discern between competing structures 
at close-packed arrangements. Varying A from 1 to 4.5, 
we have identified via analytic considerations not only 
the usual suspects for close-packed scenarios, namely hep 
and fee; also other, more complicated stacking sequences 
are observed for particular A-values - cf. Fig. [5l 

Summarizing, our investigations have given quantita- 
tive evidence about the rich wealth of self-assembly see- 
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FIG. 5: Minimum number of corona overlaps per particle 
(f7o P t) an d corresponding MECs (as labeled) for close-packed 
particle arrangements as function of A, as obtained via ana- 
lytical considerations and the MF-theory (as labeled). The 
latter gives U*(g* = y/2) = (4ttA 3 /3V2) - (1/2). 



FIG. 3: (color online). Selected MECs for the square-shoulder 
system with A = 4.5 [from (a) to (d)] and with A = 10 [from 
(e) and (f)]. MECS (a) to (d) correspond to pressure values 
indicated by arrows in Fig. [2] 
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narios of soft particles that interact via a simple, radially 
symmetric pair potential: while at high pressure com- 
pact structures are dominant, we observe at intermediate 
and low pressure low-symmetry structures, that include 
clusters, columns, and lamellae. The simple functional 
form allows to understand fully analytically the system's 
strategies to self-organize in such complex scenarios. Our 
findings can also be of technological relevance: since our 
GA-based search strategy is reliable and efficient, it can 
easily be applied to systems with more complex interac- 
tions, pointing thus towards more technological applica- 
tions, such as nano-lithography or nano-electricity. 



FIG. 4: (color online). Scaled energy per particle (U* + 
1/2) /A 3 as obtained via the MF-theory (full black line) and 
for the three square-shoulder systems investigated (as labeled) 
as a function of P*/A 3 . The vertical, dotted lines indicate the 
borders of the four different regimes. 
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